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Abstract. The phase diagram of half-doped manganite systems of formula Ao.sAo. 5 Mn03 is investigated 
within a single-orbital model incorporating magnetic double-exchange and superexchange, together with 
intersite Coulomb and electron-lattice interactions. Strong Jahn-Teller and breathing mode deformations 
compete together and result in shear lattice deformations. The latter stabilize the charge-ordered CE-type 
phase, which undergo first-order transitions with temperature or magnetic field to either Ferromagnetic 
metallic or Paramagnetic insulating phases. An essential feature is the self-consistent screening of Coulomb 
and electron-phonon interactions in the ferromagnetic phase. 



PACS. 75.30.-m Intrinsic properties of magnetically ordered materials 
tance - 63.20.Kr Phonon-electron and phonon-phonon interactions 



75.30.Vn Colossal magnetoresis- 



Manganese perovskite oxides are currently the object 
of intense activity. Motivated initially by the colossal mag- 
netoresistance phenomena, more recent studies have re- 
vealed an extremely rich phase diagram originating from 
the interplay of charge, lattice, orbital and magnetic de- 
grees of freedom Hj. The general formula is Ai^A^MnOa 
where A is in general a trivalent rare earth element (La, 
Pr, Nd) and A' a divalent alcaline element (Sr, Ca). Sub- 
stitutional doping allows to explore the full phase diagram, 
from x = to a; = 1. At the extremes, LaMnOa and 
CaMnOa are antiferromagnetic insulators. The former is 
a layered antiferromagnet, which can be explained thanks 
to the large Jahn- Teller couplings of the e g electrons of 
Mn 3+ ions j| . The latter shows a Neel ordering due to an- 
tiferromagnetic exchange of ti g electrons || . With doping, 
the double-exchange phenomena originating from Hund's 
coupling between e g and t<i g electron spins can stabilize a 
metallic ferromagnetic phase 9MlJi : Coherent band mo- 
tion occurs for ferromagnetic ordering, while strong inelas- 
tic scattering takes place in the high temperature param- 
agnetic phase. Very large magnetoresistance is obtained 
when the applied magnetic field is able to align ti g spins, 
thereby favouring the metallic phase. Nevertheless, it has 
been pointed out that spin scattering alone is not sufficient 
to quantitatively explain the phenomenon. Millis et al. 
suggested that a large electron-lattice coupling is involved, 
with the formation of Jahn-Teller polarons in the insulat- 
ing phase. Such large couplings are quite expected from 
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the very large cooperative Jahn-Teller distortions exist- 
ing in LaMnOa. Those deformations indeed involve more 
than ten per cent variations of the Mn — O bond lengths 
around all Mn 3+ ions. Local deformations have been in- 
deed revealed in charge-disordered phases by X-ray and 
neutron spectroscopy, as well as optical measurements. 
They consist of Jahn-Teller deformations around Mn 3+ 
ions, and "breathing mode" deformations with shorter 
Mn — O bonds around Mn 4+ ions. The role of these de- 
formations becomes more stringent in the charge-ordered 
phases (CO) of doped manganites. These phases strongly 
compete with the ferromagnetic metallic (FM) one at suf- 
ficient doping. Besides the Coulomb interaction between 
electrons on Mn ions, electron-phonon interaction should 
play a prominent role in this phenomenon. This is exem- 
plified by the nature of charge ordering at half-doping, for 
instance in Lao.sCao.sMnOa: While Mn 3+ and Mn 4+ ions 
alternate in two directions (say, a and b), in the other 
direction (which we here define as the c-axis), one finds 
rows of Mn 3+ or Mn 4+ ions. If CO were exclusively due 
to intersite Coulomb interaction, one would on the con- 
trary expect a Wigner crystal ordering, alternated in all 
directions. This shows that cooperative lattice distortions 
are an essential ingredient to understand charge ordering 



Charge ordering at x = 0.5 is accompanied by CE- 
type antiferromagnetic order: In the ab directions, it in- 
volves ferromagnetic and antiferromagnetic zigzag chains 
crossing each other. A qualitative explanation was given a 
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long time ago by Goodenough Q , followingthe pioneering 
structural analysis of Wollan and Koehler |J : The cooper- 
ative Jahn- Teller distortions are accompanied by orbital 
ordering, and induce the magnetic structure. Moreover, 
away from half-doping, this CE structure appears as an 
elementary "brick" to build more complicated charge or- 
dering patterns such as " stripes" [[l0| . It is thus especially 
robust and calls for a detailed explanation. 

A few models have been proposed to explain CE or- 
dering, putting the emphasis cither on intersite Coulomb 
interactions piJj , magnetism and orbital ordering . 
Mizokawa et al. juj, and Yunoki and coworkers have 
underlined the prominent role of Jahn- Teller deformations. 

Let us first list and grossly estimate the various en- 
ergy scales in the system. The on-site Hubbard repulsion 
U and the atomic level difference between the e g orbitals 
of manganese and the 2p orbitals of the oxygen are of 
the order of several eT^'s, and are larger than the total 
conduction bandwidth (W ~ 3eV). The Hund coupling 
Jh is of order leV, while the intersite Coulomb repulsion 
seems not to be larger than 0.5eV. The Jahn- Teller split- 
ting in the insulating LaMn03 phase is comparable, as 
shown by spectroscopy and optical absorption measure- 
ments In terms of a local electron-phonon cou- 
pling, it is reasonable to think of energies of the order 
of 0.2 — 0.3eV, comparable to the intersite e g hopping 
integrals to ~ 0.1 — OAeV depending on the d-orbitals in- 
volved. On the other hand, the magnetic couplings (which 
in a cubic lattice give rise to critical temperatures T c be- 
tween 100-ftT and 400K) are in the range of a few meV. 
This holds as well for the superexchange (antiferromag- 
netic) couplings as, more surprisingly, for the (ferromag- 
netic) double-exchange ones. It has been shown by Zener 
j| that T® E ~ ato, e.g. is proportional to the total ki- 
netic energy of the carriers. As will be shown below, a is 
quite small and the actual values of T® E can be easily ex- 
plained with a realistic to, for instance within De Gennes's 
mean- field picture ||. This hierarchy of energy scales is 
completed by the one set by the external magnetic field 
needed to turn the FM phase into the CE (AFCO) phase: 
It ranges from a few Teslas to 20 Teslas or more. In terms 
of energy scale per atom, this is very small, of the order 
of 0.4 — AmeV. It is thus consistent with the values of the 
magnetic exchange constants, but much smaller than all 
the other scales. This points towards an important conclu- 
sion: The stringent competition between the above phases 
require that their free energies be very close, in the range 
of a few meV per atom. Owing to the much larger electron- 
phonon and Coulomb interactions, it is reasonable to sup- 
pose that they play a dominant role in stabilizing the low- 
tcmpcrature CE phase. The necessary conclusion is that 
CE and FM phases are (meta)stablc minima of the free 
energy, separated by rather high barriers. This is consis- 
tent with the fact that the phase transitions (with tem- 
perature or magnetic field) between charge ordered and 
charge disordered phases are first-order, with strong hys- 
teresis under magnetic field. Tendencies to phase separa- 
tion between FM and CO phases have been demonstrated 
in Lao.sCao.sMnOa, Pro.7Sro.3Mn03 and other composi- 



tions. One should also notice that charge ordering is al- 
ways strong when it exists. Fine tuning of the chemi- 
cal composition between CO and FM low temperature 
phases [|17| does not allow to stabilize "weak" charge or- 
dering. This points towards strong interactions (electron- 
phonon or Coulomb) in the insulating phase, while they 
are screened in the metallic phase. This feature is over- 
looked by mean-field treatments, but can be recovered 
by taking into account exchange-correlation corrections 
to the intersite Coulomb repulsion, as shown by Sheng 
and Ting |T^] . Since the lattice distortions here also come 
from Coulomb interactions (between Mn and O ions) , we 
propose here to generalize the screening idea to electron- 
phonon interactions and use for this purpose a phenomeno- 
logical approach. 

Given the complexity of the overall Hamiltonian, here 
we restrict ourselves to a single-orbital model in two di- 
mensions, which quantitatively reproduces the various phase 
diagrams and their tuning by subtle variations of the band- 
width. Our goals are i) obtaining, for realistic values of the 
parameters, FM, CE and paramagnetic phases; ii) explor- 
ing by small variations of those parameters the different 
kinds of phase diagrams, with temperature and magnetic 
field: Of the type of Lao.sSro.sMnOs (no charge order- 
ing, FM-PM transition with increasing T) ; of the type of 
Nd .5Sro.5Mn0 3 (CE-FM-PM transitions with T, CE-FM 
with H); of the type of Pr . 5 Cao.5Mn0 3 (CE-PMCO-PM 
transitions with T, CE-FM with H). iii) obtaining first- 
order transitions between CE and FM phases. 

Taking into account explicitly orbital ordering should 
not change qualitatively the results since it works in the 
same direction ||,|l3| but may lead to quantitative im- 
provement. 

1 Model and approximations 

1.1 Hamiltonian 

According to the arguments given in the introduction, we 
assume an infinite repulsion (U — oo) between electrons 
on the same lattice site, and an infinite Hund coupling 
( Jh = oo) between the localized ti g spins and the itiner- 
ant e g spins. One can therefore consider spinless electrons, 
their spin degree of freedom being unequivocally defined 
by the direction of the local t^ g spins S. Furthermore, we 
consider in this work a two-dimensional plane of the struc- 
ture, with a half-filled band made of a single e g orbital. 
The effective model Hamiltonian is then: 

H = Hbe + Hcoui + H p h + Hse + Hh (1) 

with 

Hde =-^2 tii c \ c 3 

<ij> 

H C oul = ^ V( n i - n )( n j ~ n ) 

<ij> 
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- ^2 92Q2i(ni - n) + ^2 9bQbi(ni - n) ~ L s ^ QsiQij 

i i <ij> 

HsE — ^ [J\ — J2Qs]Si ■ Sj 

<ij> 

Hh = -9MbH ^2 S « 

i 

The first term Hue represents the double exchange 
hopping of electrons on a square lattice. Here cj , e, are re- 
spectively the creation and annihilation operators for spin- 
less electrons from a single band, and tij — ieos(0jj/2) is 
the transfer integral between neighboringMn sites whose 
ionic spins Si and Sj make an angle 9ij |p[. 

The second term Hcoui describes the Coulomb repul- 
sion between nearest neighbors (n^ = c\ci and n is the av- 
erage electron density, which is equal to 1/2 in the present 
case) . 

The third term H p h is the elastic part, which includes 
the coupling of electrons to a Jahn- Teller (JT) mode Q 2 
and of holes to a "breathing" mode Qb (#2 and gb are 
the coupling strengths, Ki and Kb the spring constants). 
In the planar geometry considered here, the other Jahn- 
Teller mode Qi is not relevant. We have also introduced 
a shear mode Q s , which is driven by Qi. Such a shear de- 
formation, which is experimentally observed at low tem- 
peratures, is essential to reconcile the alternating Mn 4+ 
breathing and Mn 3+ JT distortions wich develop in the 
ordered phases. A substantial shear deformation is in- 
deed observed in Lao.sCao.sMnOs |l9| . It results in some 
Mn — O — Mn bonds being shorter and other larger (" zig- 
zag" chains, see Fig. 1). 

The term Hse represents the antiferromagnetic (AF) 
superexchange interaction J\ between the ionic spins on 
neighboring sites, which are treated as classical. The ad- 
ditional term JiQs is a phenomenological implementation 
of the Goodenough rule: It can either enhance or reduce 
the AF coupling depending on the sign of the shear defor- 
mation, which accounts for the fact that longer (shorter) 
Mn-Mn bonds have a more (less) antiferromagnetic char- 
acter H . The last term Hh takes into account the external 
magnetic field. 

We shall study the Hamiltonian ([!]) in the mean-field 
approximation, describing the charge ordered (CO) phase 
as a charge density wave (CDW) with momentum (tt, tt). 
Let us call n A and n B the average electron densities in 
the two resulting sublattices, which correspond respec- 
tively to the Mn 3+ and Mn 4+ ions. We shall further as- 
sume that the JT coupling is only active on A sites, while 
the breathing deformations arise on B sites. With these 
approximations, the terms in the Hamiltonian which de- 
pend explicitely on (rij — n) reduce to 

H M f = -A ( n i ~ n ^ + const ^ 

i£AorB 

where the order parameter A is defined as 

A = 2V(n A - n B ) + {g b Q b + 5 2 Q 2 )/2 (3) 




Fig. 1. The CE structure and the unit cell (ions 1 to 8) in 
the ab plane: Mn i+ Jahn- Teller ions (black circles) are at the 
center of distorted Mn — O rhombus, and Mn 4+ ions (grey 
circles) undergo shear displacements, forming zig-zag chains 
with short (antiferromagnetic) and long (ferromagnetic) bonds. 

and the chemical potential has been set to zero by adding 
a term Aji = — (52Q2 — 9bQb)/2 to recover particle-hole 
symmetry (with these notations, the choice n A > n B cor- 
responds to the Q's being all positive). 

The magnetic part is also treated in mean-field, ac- 
cording to de Gennes' procedure 0], using a gaussian dis- 
tribution for the angle of the classical spins with respect to 
the mean field direction. We consider the following mag- 
netic phases: Ferromagnetic (F), paramagnetic (P), Ncel 
anti- ferromagnetic (NAF), and CE-type ordering, with 
ferromagnetic zig-zag chains, coupled anti-ferromagneti- 
cally (CE). The most general unit cell which allows to 
describe all these phases is is made of 8 nonequivalent Mn 
sites in a plane (Fig. 1). In each of these magnetic config- 
urations, the total free energy is minimized with respect 
to the following parameters: i) the magnetization on non- 
equivalent magnetic sites, ii) the average electron density 
n A on sublattice A (n B being just 1 — n A ) and iii) the 
lattice displacements. 



1.2 Phenomenological treatment of screening effects 

As mentioned in the Introduction, to have a realistic de- 
scription of the phase diagram which includes both the 
metallic and charge ordered phases at half-filling, it is 
necessary to go beyond the simple mean-field approach 
described in the preceding section by including the effects 
of exchange and correlation. Such effects on the intersite 
Coulomb repulsion in the half-doping manganites have re- 
cently been analyzed within an RPA-like calculation |Q , 
which is known to be appropriate for interacting electron 
systems at metallic densities. Since a detailed study goes 
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beyond the scope of this paper, we propose here a semi- 
phcnomcnological treatment of screening which allows a 
qualitative description of the transition between CO and 
metallic states, and which correctly reproduces the results 
of reference [|l8|. The method is further generalized to 
describe the screening of the electron-lattice interactions. 
Actually, the latter are due to the Coulomb repulsion be- 
tween Mn and O ions and are therefore also screened in 
the metallic phase. This screening should be weaker than 
that of Mn — Mn interactions, since it involves Mn — O 
rather than Mn — Mn charge fluctuations, but it should 
be sizeable. 

The procedure will be carried out in two successive 
steps. The first step consists in writing a reasonable es- 
timate for the exchange-correlation energy E xc , which is 
defined as the correction to the ground-state energy be- 
yond the Hartree mean-field result. In the second step, we 
shall define an effective hamiltonian H xc to be treated in 
the mean-field approximation, such that 



(4) 



This results in a modification of the atomic energy levels 
±A (i.e. in a reduction of the CDW gap), and it yields a 
correction to the free energy which is precisely of the form 

TP, 



1.2.1 Exchange-correlation energy 

Let us start by analyzing the simple ferromagnetic case at 
T = 0, where the electron hopping is not renormalized by 
the mechanism of double exchange. In the metallic phase, 
which corresponds to a vanishing order parameter A, the 
leading correction comes from the exchange (Fock) terms. 
These terms are responsible for an increase of the carrier 
itinerancy, which can be viewed as a renormalization of the 
hopping parameter t — » t + V(c\ci + s). Hence the kinetic 
energy is lowered by a quantity proportional to the in- 
teraction potential V, and one can write E xc = —aV (the 
parameter a is related to the dielectric constant of the sys- 
tem). On the other hand, in the CO phase, i.e. at strong 
A, the correlation energy corresponds to the interaction 
between density fluctuations on neighboring sites, each of 
them being proportional to 5n ~ t/A. Therefore, in this 
case the appropriate limiting formula is E xc ~ — V(t/A) 2 . 
These results can be generalized to the screening of the 
electron- lattice interactions, by replacing V — > g 2 /K and 
by introducing the corresponding order parameter A as 
given by eq. (||). 

A smooth interpolation between the weak and strong 
coupling behaviors can be obtained by writing the follow- 
ing expression for the exchange-correlation energy: 



E ',. 



aV + b(gl/K b + gl/K 2 ) 
l + c(f) 2 



(5) 



where a, b and c are phenomenological parameters (the 
ratio a/c — 1.44 can be deduced from ref. jy| and b/a — 
1/10 is chosen according to the ionic distances). 



As was stated at the beginning of this section, this 
formula is only appropriate in the ferromagnetic case. It 
does not account for the fact that the mobility of the car- 
riers taking part in the screening process is affected by the 
magnetic structure through the DE mechanism. We shall 
give here the arguments which allow a generalization of 
eq. (|^) to the different kinds of magnetic orderings. 

In the free-electron limit (A — > 0), where screening is 
due to coherent band motion, one expects the correlation 
energy to be reduced by a factor t/t, where 



i = i(cos(0 y /2)) 



(6) 



is the effective DE hopping parameter averaged in all space 
directions (this gives respectively 1,8/9,1/2 in the F, P, 
CE phases). The situation is slightly more complicated 
in the charge ordered phases, because E xc comes from 
incoherent hopping of the carriers to neighboring sites. 
According to Hund's rule, such processes will be allowed 
only between sites with parallel spins, which defines an 
effective number of neighbors z < z. 

In the CE phase, for instance, the lattice can be divided 
into U (up) and D (down) sites, according to the spin di- 
rection. Since each site has 2 U and 2 D neighbors, a given 
electron can only hop to the 2 neighbors with the same 
spin direction, and consequently z/z = 1/2. At finite tem- 
peratures, however, thermal fluctuations will reduce the 
absolute value of the local magnetization m. Accordingly, 
there will be a finite probability that a given U site has 
a I spin, which is given by n\j = (1 — mu)/2 (an equiv- 
alent expression holds for D sites). The total probability 
for hopping away from a U site is therefore 

2n u n u + 2n} T n} t + 2ni r n' n + 2n\ T n\ 



l u' l u 



l U ll D 



where obviously nit 



l — n\j. By adding the contributions 



for hopping processes starting from both U and D sites and 
dividing by 2, we obtain 



z = - [4 + (mu + m D ) 



(7) 



which correctly gives z/z = 1,1/2,1/2 for the F,P,CE 
phases at T — 0. Here the factors (0) and (^) introduce a 
feedback on the itinerancy of the electrons in the case of 
an applied magnetic field, which tends to align the spins 
ferromagnetically. This effect is essential in reducing the 
critical H for the CE-FM transition at low temperatures, 
as we shall see below. 

For each given magnetic configuration, instead of eq. 
(||), we shall use the following formula for the exchange- 
correlation energy: 



E xc 



aV + b(g 2 b /K b + g 2 2 /K 2 ) 



1 + ^(4) 



(8) 



where the screening parameters a, b, c have been modified 
according to 



t 

I.-. 

t 



t 



zt 
c— - 

zt 



(9) 



S. Fratini el at,: Jahn- Teller, charge and magnetic ordering in half-doped manganese oxides 



5 



We emphasize here that the terms in the numerator of Eq. 
(|J) are rescaled by the i/t factor since they arise from the 
coherent screening processes (mostly active when A — > 
0). On the other hand the local (incoherent) screening 
processes related to the term in the denominator of Eq. 
(||) are also rescaled by the effective number of accessible 
nearest neighbor sites. 



1.2.2 Mean-field potential from exchange and correlation 

We wish to define an effective hamiltonian H xc to be 
treated in the mean-field approximation such that the cor- 
rection to the free energy is equal to E xc . To this purpose, 
we replace A by an operator A (e.g. the mean-field pa- 
rameter n A is replaced by ^2i e A n t)> an d linearize the 
resulting expression. This gives 



i£AorB 



where we have defined 

B = [aV + b{g 2 /K b + g 2 2 /K 2 ) , 
The constant part in eq. (fLOf) is 



const 



(10) 



1 ' ' 7 



(11) 



- B jl + [6V(n A - n B ) + {g b Q b + g 2 Q 2 )/2] J (12) 

It is easy to verify that eq. M) holds when n A = (n A ) and 
n B = (n B ). 

One notices that a dielectric constant can be deduced 
from the screening of the gap, by writing 



VA 

A - A eff = A- Bc-p- 



which gives 



A 



1 



_ y aV+b(gl/K b +gl/K 2 ) 



2 Results 



2.1 The phase diagram: Existence of a CE phase 

The Hamiltonian in Eq. (Q) is formed by several compet- 
ing terms, and the corresponding phase diagram contains 
several phases, each one dominating in some region of the 
parameter space. To make the analysis simpler, we choose 
to vary together those parameters having similar physical 
effects. In particular, the electron-phonon couplings gener- 
ically reduce the electron mobility and, at mean-field level, 
tend to give rise to a staggered charge ordering, acting sim- 
ilarly to the n.n. electron-electron repulsion V . Therefore 
in varying V we keep constant the ratio V/(g 2 /K). For 



the sake of simplicity, we also keep a fixed J%/ J\ ratio, 
although it is not the only possible choice. 

Figure 2 reports typical phase diagrams at various tem- 
peratures as a function of the magnetic coupling J = 
J\S 2 /t with Ji — J2 (S — 3/2) and of the repulsive e-e in- 
teraction V = 0.5(g 2 / K). At low temperature (left panel) 
there is a metallic (i.e. without charge-ordering) ferromag- 
netic phase (FM) when the charge-ordering (CO) terms V 
and g 2 /K are not too large. This FM phase is naturally 
suppressed by the increase of the antiferromagnetic (AF) 
superexchange coupling J\. When the charge mobility is 
suppressed by the CO terms, one finds two distinct possi- 
ble phases. At low values of the AF coupling the pure CO 
effects dominate and a ferromagnetic (F) CO phase occurs 
at sufficiently large values of V (F-CO). The transition 
is first order, as found in Ref Q, due to the exchange- 
correlation terms. On the other hand, by increasing the 
AF coupling, the CO ferromagnetism is destabilized and 
a CE phase takes place. This latter phase naturally real- 
izes the best compromise between the electron mobility, 
favored by the ferromagnetic bonds, the CO, and the AF 
interactions increasing with J\. The CE ordering arises 
due to competing lattice displacements. In particular a 
substantial shear mode is induced in the lattice to rec- 
oncile the (asymmetrical) JT deformations occurring in 
the Mn 3+ ions with the (centrosymmetric) breathing de- 
formations around the Mn 4+ ions. The resulting lattice 
structure displays zig-zag chains formed by long bonds in- 
terlaced with zig-zag chains of short bonds. Then the pe- 
culiar CE magnetic structure naturally appears. In partic- 
ular, according to Goodenough Q , orbital ordering makes 
the sign of the magnetic couplings to be correlated to the 
length of the bonds, with AF (F) magnetic couplings cor- 
responding to short (long) bonds. Therefore, the lattice- 
driven chains with short bonds and with long bonds natu- 
rally translate into a lattice-driven CE magnetic structure. 

The temperature evolution is represented in the three 
panels of figure 2. By increasing T the weak ferromag- 
netism surviving in the CO phase at small J is rapidly 
suppressed in favor of a a CO paramagnetic phase. Also 
the CE region, being due to a delicate balance between 
CO, FM, and AF interactions, is reduced rather rapidly. 
The FM phase at small values of Ji is instead based on the 
double-exchange mechanism, which is more robust and, 
upon increasing T, is only weakly "invaded" by the CO 
paramagnetic phase. One also observes that compounds 
with a low-temperature CE magnetic structure, but laying 
near the CE-FM phase boundary, can undergo a first-order 
CE-to-FM transition upon increasing the temperature. 

It is worth pointing out that, within our formal scheme, 
the CO paramagnetic phase is the only possible mean-field 
description of an insulating non-magnetic phase at moder- 
ate temperature values. The explored temperature range 
is indeed much too low to allow for a thermal disruption 
of the CO, which instead occurs at temperatures of the or- 
der of V. However, in strong coupling, one may speculate 
that a more refined description would allow for the disor- 
dering of the charge (possibly without spoiling the local 
lattice deformations around the charges) thus producing 
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Fig. 2. Phase diagrams as function of the AF coupling Ji and the electron-phonon or Coulomb parameter V, for different 
temperatures (see text). 



the disordered paramagnetic (and polaronic) phase which 
is observed in all manganites above a few hundreds of 
kelvins. 



2.2 Sensitive role of the electronic bandwidth 

A key issue is the role of the kinetic energy in the competi- 
tion between the different phases. An extended experimen- 
tal analysis of the phase diagrams of the various mangan- 
ites [ pOpTI ] suggests that the electronic bandwidth, among 
other parameters such as lattice disorder, plays a primary 
role in determining the stability and the competition be- 
tween the FM and the insulating phases. In our model 
we investigated this relevant point by varying the bare 
hopping amplitude in front of the double-exchange term 
in the Hamiltonian (Q). We also assume that the same 
mechanism inducing the variation of the nearest neigh- 
bour hopping of the itinerant electrons in the e g orbitals 
is responsible for variations of the hopping of the ti g elec- 
trons as well. This affects the superexchange couplings J\ 
and J2, in particular Jj is expected to arise from second- 
order hopping processes of the t 2g electrons J± oc t 2 /U t2g 
(where Ut 2g is some effective repulsion between electrons 
in the same doubly occupied t2 g orbital). According to 
this assumption, when the hopping t is increased without 
changing the intersite repulsion V one moves downwards 
in the phase diagrams of Fig. 2, where the variable V/t 
is reported on the y-axis. At the same time, however, the 
increasing t produces an increase in the x-axis variable 
J cx t. Therefore, by keeping all the Hamiltonian param- 
eters fixed, but t and J\ = J2, one moves in the phase 
diagram along the dotted curves V/t — A/ J reported in 
Fig. 2. These curves correspond to similar physical sys- 
tems, where the only nearest neighbour electronic hopping 
amplitudes have been varied. 

It immediately appears that systems which slightly dif- 
fer by the electronic hopping amplitude can have differ- 
ent magnetic structures at low temperature (Fig. 2, left 



panel). In particular, by increasing t a first-order tran- 
sition can occur at low temperature from a CE to a FM 
phase. Furthermore, the temperature evolution of systems 
with different bandwidth but laying near a phase bound- 
ary can be different. This is made apparent in Fig. 3, 
where we report the temperature-magnetic field phase dia- 
grams for three systems with different hopping (and there- 
fore also different magnetic couplings J\) laying on the 
same dotted curve (V/t) = 0.04/ J. The three phase dia- 
grams correspond to systems with slightly different mag- 
netic couplings (differing at most by ten percent). Never- 
theless, already at zero magnetic field, the three systems 
display a completely different evolution in temperature. 
The more insulating (i.e., with smaller t) system having 
J = 0.037 (so that (V/t) — 1.081) never becomes metal- 
lic at zero field, but only undergoes a first-order transition 
from a low-temperature CE phase to a paramagnetic insu- 
lating phase at T ~ 0.032t (left panel, see also Figs. 4, 5). In 
the more metallic system (center panel) having J = 0.038 
(and V/t — 1.053), the CE phase disappears at a lower 
temperature T ~ 0.02i and it is replaced by an intermedi- 
ate FM phase. The ferromagnetic order and the metallic- 
ity is then destroyed at a higher temperature T ~ 0.032i 
where a paramagnetic CO phase takes place. We reiter- 
ate here that this latter phase is better to be intended 
as the mean-field representation of a disordered param- 
agnetic insulating phase. Finally, at even larger values of 
J = 0.04, corresponding to V/t = 1 the metallic phase is 
present already at low temperature and it survives up to 
aT~0.05i. 



The relevant role of the kinetic energy in stabilizing 
the uniform metallic FM phase at the expenses of the CO 
phases and particularly of the one with CE magnetic order 
is made even more apparent in the presence of a magnetic 
field. This is particularly clear in the first panel of Fig. 3, 
where the FM phase, which would be absent at zero field 
becomes the most stable solution at large enough H . 
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PrCaMnO LaCaMnO LaSrMnO 




Fig. 3. Phase diagrams as function of temperature and magnetic field, for three typical bandwidth parameters (see text). 



It is also worth mentioning that, due to the presence 
of screening, the metallic uniform solution is always a (lo- 
cal) minimum of the free energy. Therefore an (at least 
metastable) metallic solution is present even at zero field. 
The existence of a (local) minimum is a necessary condi- 
tion for the occurrence of an hysteretic behavior at the 
transition. Of course the region of the hysteresis also de- 
pends on the height of the free energy barrier between the 
maxima, of domain walls and of other non-equilibrium 
properties. Nevertheless the region in T and H where two 
minima exist provides an (excess) estimate for the hystere- 
sis region experimetally observed in half-doped Pro.sSro.sMnOs, 
Ndo.sSro.sMnOa, (Ndi-ySm^o.sSro.sMnOg and Pr .5Ca .5MnO 3 
©0- 



2.3 Discussion 

The previous subsection illustrated the main results of the 
present work: 

i) The CE phase does arise in the present one-orbital 
model and is crucially related to competition between 
the JT and breathing distortions involved in the charge- 
ordered state on Mn 3+ and Mn A+ sites respectively. The 
shear lattice deformation results from this competition 
and couples to the magnetic degree of freedom, merely 
through orbital ordering. 

ii) Exchange-correlation corrections are essential to sta- 
bilize a metallic ferromagnetic phase, due to substantial 
screening of both Coulomb and electron-lattice interac- 
tions. 

iii) The kinetic energy is a most effective parameter 
in determining the relative stability of the various phases 
upon varying the temperature and the magnetic field. 

As far as point i) is concerned, in the present work 
we show that the CE phase not only arises from elec- 
tronic mechanisms based on the presence of (at least) 
two orbitals per Mn site. The existence of a CE phase 
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case PrCaMnO of figure 3. In the CE phase, the magnetization 
is staggered. 



in a model without orbital degrees of freedom is quite re- 
markable. It is indeed repeatedly stated in the literature 
[p^|Jl^ , |2l|| that the CE phase is stabilized by the kinetic 
energy gain arising from the orbital ordering forming fer- 
romagnetic chains. Our results are not in contrast with 
this viewpoint, but underline that the above purely elec- 
tronic mechanism is not the only possible one, and that 
the coupling with lattice degrees of freedom can be of pri- 
mary importance. In this regard our results are related to 
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previous Hartree-Fock [Q and quantum Monte-Carlo Q 
calculations, where the JT deformations were claimed to 
be relevant for the occurrence of a CE phase. Our low- 
temperature phase diagrams are qualitatively similar to 
the one reported in Fig. 2(c) in Ref.Q once the distinction 
between order and disorder in the orbital degrees of free- 
dom is discarded |22]| . Our contribution in this framework 
is to show that the lattice shear deformation is a relevant 
ingredient in its own even in the absence of cooperative 
mechanisms due to electronic or JT-induced orbital order- 
ing. 

Regarding point iii) we notice that a systematic anal- 
ysis of the role of the hopping is relevant for the general 
understanding of the manganites. In the real materials 
of the general form Ri-^A^MnOs (where R and A are 
trivalent rare earth and divalent alkaline earth ions re- 
spectively) the bandwidth can be varied by changing the 
radius of the perovskite A site (where the R and A ions are 
located). Depending on the averaged ionic radius the bond 
angle of Mn-O-Mn deviates from 180° in the orthorombic 
lattice. The smaller the radius of the A site is, the larger is 
this angle, which reduces the Mn-0 overlap and the effec- 
tive Mn-Mn hopping. A systematic experimental analysis 
of this effect is reported in Ref. |jl7). The results summa- 
rized in Fig. 3 allow for a unified qualitative description of 
different half-doped materials. In particular the most insu- 
lating behavior in Fig. 3(a) is consistent with the generic 
features of Pr . 5 Cao.5Mn03. On the other hand, the cen- 
ter panel of Fig. 3 shows the same qualitative behavior 
of La .5Cao.5Mn0 3 or (Ndi_ a Sm y ) .5Sr .5MnO3. Finally 



the most metallic system in the right panel is a good qual- 
itative description of Lao.sSro.sMnOa. Nevertheless, it has 
been pointed out |p3f that a rapid change in lattice con- 
stant K, rather than necessarily small changes of t, could 
be the clue for the very different behaviours of the systems 
(A,A') .5MnO3. In our case this would correspond to an 
abrupt change along a vertical line in Figs. 2, and would 
enhance the first-order character of the transitions. 

A semi-quantitative agreement can even be reached. 
In fact, the pure double-exchange ferromagnetic critical 
temperature (J/t = 0) is, from our 2D mean-field calcu- 
lation, T c ~ 0.085<. A 3D estimate enhances this value by 
a factor 3/2 owing to the number of nearest neighbours, 
yielding T EE ~ 0.13t. For an average value t = 0.3eT^ 
one gets a transition temperature ~ 450K . It is reduced 
by the presence of the antifcrromagnetic coupling, for in- 
stance in panel c of Figure 3, in zero field T c ~ 0.05t thus 
~ 270-ftf in 3D, thus supporting De Gennes' simple mean- 
field picture. Then one obtains in the center panel the 
value T C CB ~ 18QK, and in the left panel T? E ~ 170K. 
These values are reasonable, as compared with experimen- 
tal ones, in particular one notices that T^ E is strongly 
reduced compared to T® E . This is due to the competition 
between the two order parameters (ferro and antiferro). 
Another way to understand it is to notice that in the CE 
phase the effective dimensionality is reduced by chain for- 
mation, together with charge localization this reduces the 
effective stength of double-exchange. On the other hand 
the effective antiferromagnetic exchange is close to that of 
stoechiometric CaMn03 with a T c of 12QK. 

We have also systematically investigated the role of 
the magnetic field in stabilizing the uniform FM phase. 
The typical energy differences involved in the first-order 
transitions are so small that accessible magnetic fields are 
sufficient to drive the transition from the insulating to the 
FM phases. Specifically, by taking a typical value of t ~ 
0.3eV, one can see that H/t = 0.015 (where H = gn^SH) 
roughly corresponds to ten Teslas. This value agrees well 
with the typical values experimentally used to investigate 
the (T,H) dependence of the low-temperature CE insulat- 
ing phase and the intermediate-temperature uniform FM 
phase HQ. 

3 Conclusion 

Let us compare our approach with other models which 
have been proposed to describe the half-doping compounds. 
In our treatment, the main ingredient which is responsi- 
ble for the CE-type magnetic ordering is the appearance 
of a shear deformation, with the consequent modification 
of the magnetic coupling along certain directions. In ref. 
||, two different orbitals are retained for the e g electrons, 
but the shear deformation is absent. In their approach, the 
CE order arises because the orbitals prefer to have large 
overlaps along certain directions, thus favoring the kinetic 
energy along zig-zag chains where the spins are aligned fer- 
romagnetically. In ref. |l3| , the electron-lattice interaction 
is absent, and it is again the anisotropic e g transfer ampli- 
tude of the two-orbital model which drives the CE state. 
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In both cases, however, the AF coupling J ~ O.li nec- 
essary to achieve the CE state is one order of magnitude 
higher than what is estimated from experiments, signaling 
that there must be some additional effect contributing to 
the CE ordering. Finally, we reiterate that self-consistent 
screening is necessary to explain that phases with marked 
charge order come in very close competition with metallic 
phases. We believe that this is a crucial feature of doped 
manganites, that further models addressing coexistence 
and texturing of those phases at small scales must take 
into account. 



major difference is responsible for the much lower values of 
the AF coupling Ji in our results in comparison with the 
rather large values of Jaf in Ref. ||. 
23. T. Egami and Despina Louca, J. of Supercond. 13, 247 
(2000). 
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